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Abstract. We study a many-body ground state of graphene in perpendicular magnetic fields. 
Chiral symmetry in graphene enables us to determine the many-body ground state, which turns 
out to be a doubly degenerate chiral condensate for the half-filled (undoped) case. In the ground 
state a prominent charge accumulation emerges along zigzag edges. We also show that gapless 
excitations are absent despite the presence of the robust edge modes, which is consistent with 
the Chern number C = 0. 



1. Introduction 

Graphene is a two-dimensional zero-gap semiconductor composed of carbon atoms forming a 
honeycomb lattice. Since its discovery [1], there is an upsurge of researches on graphene's 
intriguing physics, which is in may ways intimately linked to a fundamental nature in the lattice 
structure, i.e., the chiral symmetry. For one thing, the chiral symmetry guarantees a topological 
stability of the doubled Dirac cones against finite perturbations [2J. We can view this as a two- 
dimensional analogue of the Nielsen-Ninomiya theorem in the four-dimensional lattice gauge 
theory. In a magnetic field, the n = Landau level is unique in terms of the chiral symmetry, 
where the level remains anomalously sharp even in the presence of disorder that preserves the 
chiral symmetry such as ripples [3J. 

Here we ask a question: How does the chiral symmetry influence many-body physics of 
graphene? There is a prevailing view that the gap opening in the n = Landau level observed at 
high magnetic fields [4J is due to the electron-electron interaction. If we take account of the close 
relation between the chiral symmetry and the n = Landau level, the chiral symmetry may be 
again playing a crucial role in the many-body gap opening. In this paper we study the effect 
of the electron-electron interaction in undoped graphene in magnetic fields from a viewpoint of 
the chiral symmetry. 
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Figure 1. Local density of states near £^ = in the ground state 
of a half-filled, zigzag-ribbon graphene. The system consists of 
12 X 12 unit cells, and the magnetic flux is set to (/> = 1/6. A{B) 
sublattice is represented with solid (open) circles, where the size of 
each circle represents a relative value of the local density of states. 
Contribution from n < Landau levels is neglected. A charge 
accumulation on A(B) sites is seen along the top (bottom) zigzag 
edge, while the bulk has a uniform charge distribution over A and 



bottom zigzag edge B. 

2. Doubly-degenerate chiral condensate as a many-body ground state 

To model the interacting electrons on a graphene honeycomb lattice, let us adopt the extended 
Hubbard model for spinless fermions, 

H = i/kin + ^^int = t ^{e'^'^clcj + H.c) + VY^ n,nj {n, = cjc,), (1) 

where electronic spins, assumed to be completely polarized due to the Zeeman effect, are 
suppressed (hence the on-site interaction is absent). Here the screening effect is assumed to 
be so strong that we can only take the nearest-neighbor repulsion (V > 0) as the electron- 
electron interaction. The kinetic energy part is just the usual hopping, t, between adjacent sites 
(ij), and the effect of magnetic field is taken care of by the Pierls phase 9ij. 

Since it is difficult to treat the interaction exactly, we consider a projection onto the n = 
Landau level to reduce the Hilbert space, assuming that the interaction is small compared with 
the Landau level spacing. Note that in the projected subspace one can neglect the kinetic energy. 
In the presence of edges, the Landau level projection becomes nontrivial, because quantum Hall 
edge states appear between Landau levels. We identify the n = states by checking that the 
value of charge density away from the edges coincides with that for a system without edges. 

To determine the many-body ground state, we can first introduce the chiral operator F 
defined as Fc^F"^ = +c^,Fc^F~^ = —cb with F^ = 1, where c^(^) annihilates an electron 
at a site in the A(B) sublattice. H]^in has a chiral symmetry, i.e., {i^kin^T} = 0, and, for the 
one-body, £" = eigenstates tpE=o one can construct a set of eigenstates of F as (1 ± T)i^E=o- 
The even- (odd-)chirality eigenstates are localized on A(B) sublattice. If A and B sublattices 
have the same number of sites, there are the same number of even- and odd-chirality states. 
Thus, in the subspace projected to the n = Landau level at = 0, the many-body ground 
state at half-filling should be a chiral condensate, which is a state composed totally of n = 
states polarized to even or odd chirality. This is because i^int should in general have a positive 
eigenenergy E[^t > for a repulsive interaction (1/ > 0), whereas a chiral condensate exactly 
has E'int = for the nearest-neighbor repulsion. We have numerically confirmed that this is 
indeed the case, even in the presence of zigzag edges, through an exact diagonalization of i^int 
for finite systems. We stress that the even and odd chiral condensates, |^^+) and |^^-), constitute 
a doublet |^^) = (|^+), 1^-))? which implies that the ground state is not a simple charge density 
wave (CDW). The detail of the doubly degenerate ground state will be discussed elsewhere [5]. 

3. Effect of zigzag edges 

Now that the many-body ground state is determined, we can calculate the physical quantities. 
Here we consider a zigzag-edges ribbon with a periodic boundary condition along the ribbon, 
and examine the influence of edges on the ground state. To see this, we first show the charge 
distribution in the ground state in a zigzag ribbon. In the ground state with completely filled 




Figure 2. (Color online) (a) For a zigzag ribbon, schematically shown in (a), the charge-charge 
correlation function along the zigzag edge direction is displayed in (b) as a function of distance 
between two sites for several distances from the edge; (c), 1 (d), 2 (e) and 3 (f). The sample 
has 36 X 36 unit cells with half the sample size of 36 (with each unit cell containing two sites), 
and the magnetic flux is set to (j) — 1/12. The correlation function is normalized by square of 
the bulk charge density (n) = 1/2. 

n < levels, charge density shows no symmetry breaking with a uniform value 1/2. If we focus 
on the n = level, however, the charge distribution exhibits characteristic patterns near the 
edges as shown in Figure [TJ which should be observed experimentally as a local density of states 
(LDOS) around E — {) with a scanning tunneling microscope. Here the size of a circle denotes 
the relative value of LDOS for the n = level. The magnetic flux penetrating through a unit 
hexagon is set to = 1/6. One can see that, near each edge, there is a charge accumulation, 
which is localized on A (B) sublattice for the top (bottom) edge. In this sense a CDW (or 
sublattice selected charge accumulation) is induced by the edge. This behavior reminds us of 
the presence of one-body edge modes with a flat band near a zigzag edge [6J. In the bulk region, 
on the other hand, the CDW decays and the value of LDOS converges to the uniform bulk value. 
It is easy to check that the CDW decays on the magnetic length scale Ib away from the edge, 
where Ib = 3^/^a/\/27r0 ^ 2a for the present value of with the lattice spacing a. 

Next we examine the charge-charge correlation along the edge direction. In Figure [21 the 
two-point correlation function is plotted as a function of the distance between sites i and j. 
The correlation function is calculated along zigzag lines enclosed with dashed loops depicted in 
Fig. 2 (a), and the panels (b) show how the charge-charge correlation changes as we go away from 
the edge. Sample size dependence is hardly observed between 24 x 24 and 36 x 36 unit cells. 
Directly on the edge (c), the charge-charge correlation oscillates with a period of two, indicating 
realization of the A-site selective charging with a long-range order. The large amplitude of the 
oscillation reflects sharp accumulation of charge right on the edge. If we go away from the edge 
into the bulk(d-f), the oscillation diminishes, but the charge-charge correlation is preserved, with 
a saturation over a distance from the edge as small as 3. Thus the many-body ground state is 
shown to exhibit a CDW-like behavior over the whole system. We should stress, however, that 
the CDW's in the edge and bulk regions are essentially different, since the former is a consequence 
of the zigzag edge state locahzed on one sublattice, while in the latter the doubly-degenerate 
ground state still preserves sublattice symmetry. 

This behavior is reminiscent of the v — 1/q fractional quantum Hall state in two-dimensional 
electron systems with an odd q. In the latter case, the ground state is a quantum liquid with a 
topological degeneracy Nb, — q'ln the bulk, whereas the state is locked to an edge state in the 
edge region. By contrast, the present ground state is different from the fractional quantum Hall 



liquid, in the sense that here the topological degeneracy = 2 is a direct consequence of the 
chiral symmetry, and fate of CDW is strongly influenced by the shape of edges. 

4. Chern number 

Finally we examine topological properties of the zigzag ribbon in the ground state by calculating 
the Chern number. Exact diagonalization of i^int gives a flnite gap between the ground state and 
the flrst excited states, so that the Chern number calculated for the ground state is well-deflned. 
Note that the Chern number should be formulated in a matrix (non-Abelian) formalism [5], 
since the ground state is now doubly degenerate. Thus we calculate the non-Abelian Chern 
number |7j as 

where A is the Berry connection. The many-body ground state at half-fllling turns out to 
have C = 0, which indicates that the ground state is a Hall insulator. While one may have 
a naive expectation that the presence of edge states is accompanied by gapless excitations, 
both of the energy spectrum and the Chern number suggest the absence of gapless (conductive) 
excitations. This may seem contradictory, but is actually consistent with an observation that, 
while a topological state (with a nonzero C) implies an existence of edge states, converse is not 
necessarily holds. 

5. Summary 

We have studied many-body states in the undoped graphene in a quantum Hall regime, and 
determined the many-body ground state utilizing the chiral symmetry. The ground state is a 
doubly degenerate chiral condensate, where the chiralities of the fllled states are fully polarized. 
Investigation of charge density and charge-charge correlation has revealed that the ground state 
has charge accumulations with a selected A (or B) sublattice emerging along zigzag edges. A 
CDW-like behavior is observed also in the bulk, but its origin is essentially different from the edge 
case. Exact diagonalization of i^int shows that there is no gapless excitation above the ground 
state despite the presence of the robust edge mode, which is consistent with the calculated Chern 
number of zero. 
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